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Abstract 

In trying to explain the anomalously sharp dependence of thermal conduc- 
tivity on magnetic field in the cuprate high t c superconductors, it was found 
that previous discussions of quasiparticle motion near the gap nodes were in 
error because of failure to take into account total current conservation. We 
present corrected equations of motion for quasiparticle motion and discuss 
their relationship to the observations. 
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In thermal conductivity of High T c cuprates, at all temperatures below T c , there is 
observed a remarkably steep dependence on applied c-axis magnetic fields even as low as 
.IT. Fig (1) shows observations in both (La — Sr^CuO^ and YBCO |I[ From thermal Hall 
(Righi-LeDuc) effect, it was concluded that the enhancement of thermal conductivity below 
T c is electronic, and that large fields seem to block this conductivity. But one of the primary 
mysteries has been the very sensitive dependence on magnetic fields in which any Zeeman 
or cyclotron frequency corresponds to an energy < 1°K, much lower than thermal energy 
~ 30-60°K in either case. 

It has been clear that the carriers responsible for electronic thermal conductivity at 
lower temperatures must be those near the nodes of the "d— symmetry" energy gap at ±kp, 
±fcp. Although their density is relatively small, the very rapid removal below T c , as the 
gap opens up, of the electronic scattering mechanism responsible for the normal state "T- 
linear" resistivity could allow peaking of quasiparticle conductivities. 0] The ARPES data on 
electron Green's functions near the gap nodes, however, do not support this interpretation, 
and it seems more plausible that the extra thermal conductivity is caused by the anomalous 
properties of carriers near the nodes — specifically, that they acquire velocitiy parallel to 
the Fermi surface as well as the conventional Fermi velocity. This problem will be the 
subject of later publications. Our first hypothesis for the i7-field sensitivity was that spin- 
charge separation was somehow involved in this large conductivity enhancement and that 
k was being destroyed by the effects of Zeeman splitting, but this mechanism seems unable 
energetically to account for the sharp B- dependence. The mechanism is described below 
and involves quasiparticles which are almost conventional. 

This problem has been discussed by Lee and Simon ||, who derived scaling laws for 
thermal conductivity which should be valid once the carriers are confined to the nodes. This 
should be true for T < T c /2, at the very least. These scaling laws are 




(la) 
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and 

K xy = const x£f' (^?) . (16) 

These scaling laws are quite successful for T < 50° for optimally doped YBCO, and 
clearly the arguments of Ref (3) (which depend on the result that energy levels in the nodes 
scale as y/~B, and hence depend sensitively on magnetic field) contain the essence of the 
phenomena. The purpose of this paper is to support their general scaling argument with a 
specific picture of the nature of the energy levels at the nodes in a magnetic field; and to 
predict that broadened Landau-like levels will be observed, contrary to previous discussions 
of these energy levels which ignore the subtle properties of gap nodes. 

One can discuss the physics semiclassically either as a Tomasch interference effect in- 
volving Andreev scattering or as a form of "Landau level" formation for excitations around 
the "Dirac point" nodes of the <i-wave gap. The former point of view is simpler. First we 
consider the effect of the uniform magnetic field and will later discuss the effect of vortices. 
The B field (along the c-axis, perpendicular to the plane) causes a Larmor precession of all 
the electrons around the Fermi surface according to 

g 

hk = - vf x B (electrons and holes) (2) 

Here e is the electron charge, negative in this case, vf the Fermi velocity at the nodal point, 
and hk the electron momentum. Since we speak of an electron, \k\ may be assumed to be 
somewhat (5k) greater than k F , giving an excess kinetic energy 5E = v F 5k. Holes obey 
the same equation with the signs of both e and vp changed, hence with the same direction 
of precession. The fc-vector of an electron near a node is rotated away from the node, and 
therefore eventually encounters a gap equal to its excess kinetic energy. We assume that in 
the superconducting state there are point nodes in the gap, which in this region varies as 

A = ^ (A;_l - k nodc ) = hv A (k± - k node ) = hv A Ak . (3) 

where k ± is the momentum perpendicular to v F = |§ and hence parallel to the Fermi 
surface, and v A is defined by (3). At this point the electron must be Andreev reflected as 
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a hole in the state of opposite momentum, supplying (— 2e) of charge to the condensate. 
(Essentially, Fermion number as well as charge must be conserved.) This hole Larmor 
precesses in the same direction in /c-space. 

This at first appears to scotch the possibility of a periodic process. However, it must be 
realized that the Andreev reflection takes place not on a zero-momentum condensate but on 
a condensate which macroscopically obeys London's equation 

eA 
p = — 



so that the momentum carried off by the pair is ^p. This means that the hole reappears 
not at the same momentum but at a momentum shifted back by 



Sp = -2h A k 

(it is easily verified that the change in ^f- = 2UAk.) Now the hole Larmor precesses from 
one barrier to the other, is again Andreev reflected as an electron, and the process begins 
again. 

The period of this process defines a frequency / which will be the excess energy level 
separation for electrons of that energy. It is easy to estimate. For electrons of energy hu, 
the separation of the two barriers is 

dk 



Ak 



dA 



(4) 



and the velocity in /c-space gives 



dk± 
~dt 



f 



Ak 1 evp 
~Y ~ 2~hc 



■B 



or, 



/ 



Acuj 



dA 



dk 



B 



To find the lowest energy state permitted by the uncertainty principle , we set / = |£ and 
obtain 
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n 2 ul =e^ VfBva = \n 2 {^) 

= f (H) E F ■ % 

(5) 

We have defined a velocity va = or to represent the steepness of the energy gap at the 
nodes. This has been estimated by Lee j| to be ~ \ V F- Using that figure, and E F ~ lev, 
at IT Hujq ~ 4.5mev, which is of the right order of magnitude to explain the data. 

A very physical way to think of this resonance is as a Landau level in the Dirac points 
at each of the four nodes of the energy gap, at which the quasiparticle excitations act like 
massless Fermions of energy 

El = tfv%{ h -k F y + (^j kl (6) 

These excitations change character as they precess around the elliptical energy contour, from 
electrons through neutral electron-hole mixtures to holes, but by the same argument given 
above the quasiparticle precesses everywhere as though it had charge e and a velocity given 
by 

vq p = t V fc E k , (7) 

i.e, 

h k = ^xB . (8) 

As with the Andreev reflection picture, (8) seems a necessary consequence of basic conserva- 
tion laws. So long as the Larmor frequency remains small compared with the maximum gap 
Ao ~ 50mev tunneling through the gap maxima will play a minor role. (8) leads to the same 
result, (5), for the cyclotron orbit frequency except for constants of order 1. Although in 
momentum-space the orbit is slightly complicated by the momentum shifts due to Andreev 
reflection on a moving condensate, these shifts occur only parallel to the Fermi surface, and 
as far as real space is concerned it is a conventional elliptical orbit. 

In order to go beyond the semiclassical picture, it is necessary to solve the actual 
Bogoliubov-de Gennes equations for the quasiparticle wave function near one of the fc-space 
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nodes of the gap. These are 

Hip = huo^ (9) 



where ip is a two-component vector and H a two-by-two tensor 



n = 



^ vf ■ (p — -A — mv s ) — fj, A* 



(10) 



\ A (-v F ■ (p + \A + mv s ) + n) J 

We will specialize to a single node of the gap and rotate this node to k x = k F , k y = 0. 

We have included the superfluid velocity in (10), thus reintroducing the vortices. The 
superfluid velocity is equivalent to a displacement of the entire Fermi surface (even near the 
nodes) by mv s = p s in momentum space, which leads to a "doppler" shift in the energy 
relative to the chemical potential. Its sign agrees with that of A. 

But it is also essential to include the position dependence of A in the presence of a 
magnetic field, which is not irrelevant because the off-diagonal, anomalous terms of (10) 
represent processes in which a particle Andreev scatters into a hole, with the charge of the 
missing pair going into the condensate. Such a process must also be momentum-conserving, 
so that if the condensate has been accelerated by the magnetic field the self-energy scattering 
must reflect that fact, and the pair created in this process must carry the momentum of 
condensate pairs. 

What is the space-dependence of A? We can get a reasonably good characterization in 
the "London" regime H cl « B « H c2 where the vortices are dense enough to overlap 
strongly so that the field is fairly uniform, while the vortex cores still comprise a very 
small relative volume oc Thus almost everywhere the field and current satisfy London's 
equation as modified by the presence of a vortex lattice. A is constant in magnitude, while 
its phase ip satisfies 

j s = n s ev s = — (HWtp - —A) . (11) 
m c 

j s is zero averaged over the interior of the sample, where the field is nearly uniform. But 
near a vortex j s behaves like - and ip = 9, the actual angle measured around the vortex 
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core. Note that A is single- valued everywhere except at the vortex cores, though ip of course 
is not. 

Now let us insert (p into the Hamiltonian (10). There are two ways to transform away 
the spatial dependence of A: to transform unitarily by 



U 



V 



o 



n' = u^nu 



(12a) 



or by 



U 



v 

(These are the only two single-valued transformations) In the first case, 



(126) 



n 



holes 



\ 



(13a) 



v F ■ [(p + -A) + mv s ] A 

A* Vp ■ [p + \A + mv s ) J 

using TiVtp = mv s + —A. This is the Hamiltonian appropriate for hole-like excitations, since 
when A vanishes at the node, the wave-function is purely hole-like. For electrons, we use 
(10b), and obtain 



electrons 



V 



(136) 



v F ■ (p- - c A- mv s ) A 

A -vp ■ [(p - - c A) - mv s 

If we neglect the (pp ■ v s ) term, we come to the conclusion that the basic form of the 
wave function is cyclotron orbits in momentum space around Dirac points at the gap nodes, 
with frequencies proportional to \[B. 

The fact that these are two separate equations, one for electrons and one for holes, reflects 
the fact that charge is actually exactly conserved: when an electron is Andreev scattered 
into a hole, the charge goes into the condensate and is recovered in the inverse scattering. 
As Volovik has emphasized, gap nodes are sites of chiral anomalies where supercurrent 
and quasiparticle current are not separately conserved. 

But what of the pp ■ v s terms? These are quite strong in a localized region around the 
vortex core. We predict that these will have only a moderate effect. The cyclotron orbits of 



low order have wave-functions which occupy roughly an area of n (magnetic lengths) 2 , which 
is of order the size of the vortex lattice unit cell, and it will cost energy to localize them 
further. I suspect that the vortices act as moderately strong scattering centers, broadening 
the Landau levels without changing them quantitatively (Since % — 0, these potentials will 
average to zero.) 

Throughout this discussion, we have chosen a gauge such that A y — 0, A x — By. The 
dependence of A on k y is then straightforward. But if we are to assume general gauge 
covariance, clearly 

Afo)=A(A*-^A) 

This in turn, implies that the v s term will also enter in the y direction by symmetry. Gather- 
ing all of those various strands together, we come to the wave equation which nodal particles 
(electrons, for example) obey, 

= Hi) 



n 



^ Vp ■ [p — - c A — mv s ] i>a ■ [p — ~A — mv s ) ^ 



(14) 



\ • (p - \A ~ mv s ) -v F -(p-- c A- mv sj j 
Thus the particles at nodal points are, unexpectedly, in all respects equivalent to Dirac 
fermions in a vortex lattice and a magnetic field. At this point, the equations have become 
essentially equivalent to those discussed by Volovik |J in calculating Fermionic energy levels 
around the gap nodes in the He — 3 A phase. Volovik incorrectly asserts that the equations 
for a superconductor in a magnetic field are not equivalent to his equations, for which the 
effective magnetic field, "B"= V x I, is provided by the order parameter texture; but as 
we have shown above, this is not the case near a gap node and the equations are identical 
except for the doppler term v ■ P s , as is required by conservation laws. 

In the absence of this term, the eigenvalues are easily obtained by squaring the Hamil- 
tonian. Following Volovik, we write 
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c ■ 



where f is the Nambu "isospin" matrix homologous to Pauli spin. Squaring this, we obtain 

+ ^v f v A r z B z 



(15) 



This leads to a spectrum of eigenvalues 



E 2 n = h 2 v F v A ^)[(2n + l)\e\-r z e) (16) 

The striking feature of this is the "zero mode" n = 0, r = A, which is caused, basically, by 
the topological effect of the vortex cores on the quasiparticles, and in fact there is just one 
state per vortex for holes and one for electrons. 

What do we predict, from this picture, for the thermal conductivity? The physics is that 
of electrons in Landau orbits whose spacing is proportional to \[B 1 being scattered by vortex 
cores (also spaced by Ib oc VB), either incoherently (at low fields, and high temperatures, 
where the higher orbits are relevant and the lattice is disordered) or possibly coherently, 
at high fields and low temperatures. In the latter case there will be gaps in the spectrum 
and k could fall to near zero: this is likely to be the cause of the plateaus observed in some 
samples. 

The main predictions do not differ very much from those of Lee and Simon's successful 
scaling laws, at least at low temperatures. 

As long as we can ignore the vortex cores the same conclusion they reached, that the 
entire energy level structure scales as v^B, seems entirely valid, so that the scaling function 
will be F(\fB /T). (Just by dimensional analysis: changing B is equivalent to a shift in p by 
a factor .) The justification given here is slightly more rigorous than that in the reference. 
The function which Ong and Krishana use for K xx , 

1+pB 
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obeys the Lee-Simon scaling law over most of the range where it should, p having roughly 
the form 1/T 2 ; and as to order of magnitude, pB is of order unity where the lowest cyclotron 
orbit ujq is of order T. 

At low fields and moderate temperatures, we could imagine that the classical formula 

k oc To oc T 1 (17) 

1 + ijJ l c T l 

would be usable. uj"1 is proportional to B, and it is reasonable that 1/r caused by scattering 
on the cores would be independent of B and simply proportional to the density of final 
states oc T (though by no means obvious.) At high fields the states become increasingly 
localized and eventually only the zero mode is occupied. It is not at all clear why cra-^ in 
this region and in fact plateaus with o ~ are occasionally observed. (There are actually 
more low-energy states, not fewer; but they seem likely to localize better as B increases.) 
The experiments show us that the vortex current scattering effect is, as expected, similar 
in magnitude to the cyclotron frequency itself, so that the cyclotron orbits are not sharp in 
energy except where the scattering is coherent. 

As for K xy , as pointed out by Lee and Simon this is smaller than k xx by a factor ~ T/Ep 
because it comes from hole-particle dissymmetry. It is significant that over a range of 
intermediate -B's it fits fairly well to \[Bj\ +pB, (see Fig. (2),) which would follow from 
semiclassical arguments from the Landau orbits. But since this is in principle outside the 
"node" approximation, it will not be discussed in detail here. 

The "spin gap" phase seems to be one in which the mode structure of the pseudogap 
is already developed. (This is the picture proposed by Lee and Nagaosa ||, following 
Baskaran and Zou, and borrowed by Fisher, Balents and Nayak ||.) Thus some of the 
above considerations may be relevant for critical phenomena near T c , in the spin gap phase. 

I would like to acknowledge very extensive discussions of data and theory with N.P. Ong, 
as well as the cooperation of K. Krishana in making his extensive data readily available. 
Discussions with C. Nayak, R.B. Laughlin, and F.D.M. Haldane were very stimulating. 
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